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Abstract 

A nonlinear system possessing a natural forbidden band gap can trans- 
mit energy of a signal with a frequency in the gap, as recently shown for 
the nonlinear chain of coupled pendula [Phys Rev Lett 89 (2002) 134102]. 
This process of nonlinear supratransmission, occurring at a threshold ex- 
actly predictable in many cases, is shown to have a simple experimental 
realization with a mechanical chain of pendula coupled by a coil spring. It 
is then analyzed in more detail by first going to different (non-integrable) 
systems which do sustain nonlinear supratransmission. Then a Josephson 
transmission line (one dimensional array of short Josephson junctions cou- 
pled through super- conducting wires), is shown to sustain also nonlinear 
supratransmission, though being related to a different class of boundary 
conditions, and despite the presence of damping, finiteness and discrete- 
ness. Finally the mechanism at the origin of nonlinear supratransmission 
is a nonlinear instability, and it is briefly discussed here. 



1 Introduction 

It has been recently demonstrated that, in addition to energy spectral localiza- 
tion (the Fermi-Pasta-Ulam recurrence phenomenon PQ) and to energy spatial 
localization (soliton generation from initial data [U ED a nonlinear chain 

of oscillators possess another striking fundamental property called nonlinear 
supratransmission |jy. This phenomenon was shown to occur in the nonlinear 
sine-Gordon chain, which possess a natural forbidden band gap, when it is sub- 
mitted to irradiation at a frequency in the stop gap. While in a linear chain 
the signal would exponentially vanish in the medium, it does not do so in the 
nonlinear case if its amplitude exceeds a threshold value. 

In |H] the irradiation of the medium was modelized by prescribing the bound- 
ary value at one end of the chain. It is worth mentioning that nonlinear trans- 
mission does occur also in the case of a true wave scattering, namely when a 
monochromatic plane wave scatters onto a nonlinear medium with a frequency 
in the gap 0. 

•Journal of Physics Condensed Matter 15 (2003) 2933 
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In the case of the sine-Gordon model, the threshold can be predicted exactly 
by invoking the static one- breather solution 6 . In short, energy penetrates the 
medium as soon as the amplitude A of the harmonic driving, at a frequency f2 in 
the gap, exceeds the maximum amplitude of the static breather of frequency f2. 
This energy travels then through the medium by means of nonlinear localized 
excitations (kinks, breathers, solitons). 

The purpose of this paper is, after recalling results of 0, to display experi- 
mental results on the mechanical chain of pendula and to discuss the generality 
of our result by exploring extensions to different situations. We shall in par- 
ticular explore the robustness of nonlinear supratransmission when the medium 
experiences damping, different nonlinearities and different classes of boundary 
values. 

The next section is devoted to a short reminder of the results published 
in [H] and it is intended to settle formalism and basic facts about nonlinear 
supratransmission. The model is the sine-Gordon chain submitted to a Dirichlct 
condition at the origin on a vanishing initial background (pendula at rest). 
Then section [21 relates the experiment made with a mechanical pendula chain, 
coupled by means of a coil spring, and which is forced at one end by a periodic 
torque. A systematic exploration of the chain response in a frequency range 
within the gap shows spectacular agreement with the theory. In section 0] we 
numerically describe the characteristics of the breathers generated in the sine- 
Gordon chain and discover some simple relations between the parameters of the 
emitted breather with respect to those of the boundary driving. The section 
El deals with the energy transmitted by the nonlinear medium as a function of 
the driving amplitude at given frequency in the gap. Particular emphasis will 
be put on the effectiveness of the effect in different, not necessarily integrable, 
cases. The nonlinear instability which is the generating mechanism of nonlinear 
supratransmission is briefly discussed in section Although the mathematical 
analysis is still to be constructed we propose a quite simple illustration of the 
process by a perturbative analysis of the sine-Gordon system driven close to a 
breather mode. 

As another domain of study, we consider in scction[7|the sine-Gordon model 
for Neumann conditions at the boundary (the derivative at the origin is pre- 
scribed). This is a model for a chain of short Josephson junctions whose first 
one is submitted to an external AC current. By using numerical simulations, 
the process of nonlinear supratransmission is shown to hold, and the threshold 
of energy transmission to obey a similar simple rule. Remarkably, in this case 
the energy flows by means of kink (or anti-kinks) and not by breathers (or kink 
anti-kink pairs). This is an interesting issue in view of applications as kinks 
are the objects that have an experimental signature (through the Josephson 
current). 
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2 The sine- Gordon chain 



2.1 Generalities 

The following normalized discrete sine-Gordon chain of N locally damped cou- 
pled oscillators u n (t) 

u n - c 2 (u„+i - 2u n + Un-i) + sinu„ = -j(n)u n , n=l,--- ,N , (1) 

is submitted to the boundary value at the origin (Dirichlct condition) 

u (t) = A sin fit, (2) 

acting on a medium initially at rest, namely 

u n (0) = 0, u n (0) = 0. (3) 

In this section the damping coefficient 7 is used to model a semi-inhnite chain 
by an absorbing boundary. More precisely we take ujy+i(t) = with 

> n < m : 7(71) = 0, 

, , 2n — m — N . , 
m<n<N : j(n) = a[l + tanh( — )], (4) 

where the parameter b is adjusted to have a damping factor 7(71) varying slowly 
from almost to almost 2a on the last (N — m) particles. A typical experiment 
will have e.g. N — 100, m — 60, a = 0.5 and 6 = 3. Note that in section 
we shall use instead the reflective boundary condition un+i = un for an 
homogeneous damping 7(^1) = cste for all n. 

The equation Q is considered as a second-order ordinary differential system 
for the N coupled oscillators u n (t) . This system is then solved with the subrou- 
tine dsolve of MAPLE software package which uses a Fehlberg fourth-fifth order 
Runge-Kutta method. 



2.2 Forcing in the gap 

The linear dispersion relation ui{k) of the chain is 

lu 2 = l + 2c 2 (l-cosfc). (5) 

For a driving boundary J3J) with frequency in the forbidden band gap (FBG), 
namely fi < 1, a linear chain would sustain the solution (evanescent wave) 

Asin(fi7;)exp[-An]. (6) 

The parameter A is given from the dispersion relation JSJ written for ui = fi < 1 
and k = iX, namely 

1-fi 2 ' 



A = arccosh ( 1 H — — ) . (7) 
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In the nonlinear case, in order to fit the boundary condition (J2J, the ap- 
proximate solution (exact in the continuous limit) is given instead by the static 
breather 

Ac sin(Slt) 



Ub(n,t) = 4arctan 
where the breather center no solves 

A = 4 arctan 



SI cosh(A(n — no)) 



Ac 



SI cosh(Ano) 



(8) 



(9) 



The above equation expresses that the static breather fits the driving field J2Jl 
by adjusting its position such as to match the amplitude A. The spectrum of 
the breather signal does not match exactly the monochromatic forcing J2J) and 
it adapts by sending phonons at third, fifth, etc., harmonic frequencies. 

The above assertion is checked by performing two simulations of P at a 
given forcing frequency, say SI = 0.8, and a given amplitude, say A = 2, in 
one case with the harmonic forcing (J2J) , in the other case with the breather-like 
forcing Ub(0,t) where the value of no is calculated from Q with A = 2. To 
avoid initial shock we also set initial velocities matching the ones of the static 
breather. Then we evaluate by Fourier transform the spectra of one particle 
of the chain in both cases (we have selected the particle 50 on a chain of 100 
pendula). The result is displayed on figure ^ where we see that phonons at 
frequency 3S1 are indeed emitted for harmonic forcing while no phonon appear 
for a breather-like forcing. The small contribution at frequency 1 in figure (a) 
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Figure 1: Spectra of the particle 50 submitted to, (a) harmonic boundary driving 
at frequency 0.8 and amplitude 2 and to, (b) breather boundary driving at 
frequency 0.8. 



represents a collective motion of the pendula generated by the initial shock 
resulting from inadequation of the initial velocities with the boundary value. 



2.3 Bifurcation threshold 

The adjustment of the breather center no provided by 10 not always has a 
solution. Indeed it works up to the maximum value A = A s of the breather 
amplitude realized for no = 0. Beyond this threshold, for a driving boundary 
with A > A Sl we have shown in |S] that nonlinear supratransmission occurs, i.e. 
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the medium starts to transmit energy by means of nonlinear modes generation 
(breathers and kink-antikink pairs). From (JSJ, the threshold A s reads as the 
following function of the frequency f2 



A, 



4 arctan 



(10) 



which has the approximate value 



A s ~ 4 arctan[ 



n 



(ii) 



which would hold for the continuous sine-Gordon equation obtained in the limit 
1/c — > 0. Note that the maximum difference (when O varies in [0, 1]) between 
the above two expressions is already 1% for c = 1 and goes down to 0.01% for 
c = 10, the value we have used in most of the numerical simulations. 
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Figure 2: Representation of u n (t) as a function of time for n = 60 in the case 
SI = 0.90 for two amplitudes. 



This qualitative definition of a bifurcation threshold can be checked on nu- 
merical simulations of Q with the boundary condition J3J) by varying, at given 
frequency f2, the amplitude A around the above value A s — 1.803 (for c = 4). 
There are many means to determine appearance of nonlinear supratransmission, 
a simple one being the observation of the motion of one particle of the chain. As 
an example we display in figure [3 the motion of the particle 60 of a chain of 200 
particles driven at frequency 0.9 at amplitudes A = 1.78 (no supratransmission) 
and A = 1.79 (supratransmission) for a coupling factor c 2 = 16. Each large os- 
cillation in the second figure corresponds to a breather passing by. Two of them 
are generated and cross the site 60 at times 120 and 160. The small oscillation 
seen between the humps are the harmonic phonons, mainly of frequency 3f2. 

For illustration we show on figure 0a picture of oscillators amplitude u n (t) 
at given time (here t — 120) in the case of nonlinear supratransmission obtained 
for A = 1.8 and £1 = 0.9. A first breather propagates to the right while a second 
one is just being generated near the origin. Note the amplitude of the breather 
with respect to the one of the driving boundary. 
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Figure 3: Plot of u n (t) as a function of n at t = 120. 



By a systematic exploration of the chain response we draw figure0]obtained 
for 200 particles with a coupling c 2 = 100 (some experiments have been actually 
made with smaller coupling and less number of points to shorten computation 
times) for a typical time of 200 (for frequencies close to the gap value 1, time had 
to be increased up to 500). The points on figure0]are obtained with an absolute 
precision of 10 -2 for the amplitude A. They are compared to the theoretical 
threshold expression (|10fl (continuous curve). 




Figure 4: Bifurcation diagram in the (A, fi) plane. The solid line is the threshold 
expression <|10fl ■ The crosses indicate the lowest value of A for which nonlinear 
supratransmission is seen on numerical simulations. 

The figure 01 shows excellent agreement to formula (f 101) apart for some dis- 
crepancies starting below 0.33 and 0.18. This results from the driving which, 
thanks to the nonlinearity, generates phonons at multiple frequencies (here third, 
fifth, etc.). If these frequencies lie in the phonon band, the phonons move away 
from the boundary and have no effect on the forcing. If however they lie in the 
FBG, the related phonons do not propagate (which we call phonon quenching) 
and stick on the boundary where they add contribution to the driving. This 
effect indeed disappears when driving the system Q with the exact breather 
boundary value Uf,(0,i) for which we have checked that nonlinear supratrans- 
mission never occurs at an amplitude A < A s , while it occurs already at A = A s 
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(or for very small deformation of the perfect breather). 

3 Experiments on a mechanical chain 

Our purpose here is to show that the process of nonlinear supratransmission can 
be easily realized experimentally on a chain of coupled pendula as depicted in 
figure [5J This chain has been built following M. Remoissenet [S]. The pendula 
rotate freely around a piano wire stretched between two supports, they are 
coupled together by a coil spring tightened to each pendulum by a screw. 



15 mm 7 mm 




8 mm 



Figure 5: Sketch of the mechanical chain. 

Such a chain of 48 pendula is driven by an electrical engine steered by a 
generator of sinusoidal tension. Upon varying the frequency, at low amplitude, 
in the phonon band, we can determine the parameters of the chain (coupling 
constant a 2 and angular eigenfrcqucncy Vq) by measuring the wavelength 2-k/U 
of produced wave (in units of pendula number). By comparison to the dispersion 
relation of the linear chain 

v 2 = vl + 2<r 2 (l -cosfc). (12) 

we determined the following parameter values 

i/ = 15 Hz , a = 32 Hz (13) 

Then in a time normalized to the eigenfrcqucncy, we have the dispersion relation 
JSJ) with 

c= — = 2.13, (14) 

which is the fundamental parameter of the model, entering in particular expres- 
sion (|l(Jfl of the threshold amplitude for nonlinear supratransmission. This is 
the formula we want here to confront to experiments. 

We proceed with a systematic exploration of the chain response to a signal 
frequency in the gap (of value less than vq). The method consist in increasing 
slowly, at fixed frequency, the amplitude of the driving, up to the time when a 
nonlinear mode is seen to be generated (a picture of a breather generated at a 
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Figure 6: Experimental values of nonlinear supratransmission threshold com- 
pared to expression itTUl) . 



frequency of 0.85 in normalized units is displayed in Repeating 3 times, for 
each driving frequency value, the measurements of the driving amplitude that 
generates a nonlinear mode, we eventually obtain the figure [S] which displays 
the measured threshold amplitude in terms of the signal frequency (normalized) . 
There, the full line curve is the function A s (£l) given in (|10fl with c = 2.13. 

Despite the small number of pendula, inducing reflection at the open end, the 
inherent damping, and other mechanical imperfections, the experiments provide 
a spectacular realization of the theoretical threshold prediction. 



4 Characteristics of the emitted breathers 

One important issue concerns the nature and the characteristics of the nonlinear 
modes which propagate in the medium. In the first instance, we have to check 
that the emitted structures correctly match the moving breather, at velocity 
v < 1, which is convenient to write as 



u v (n,t) =4arctan 



' / c 



Vl+r 2 Vl-n 2 



cosh 



if ' c — vt) 



Vl+r' 2 Vl-v' 2 



(15) 



By choosing correctly the two parameters r and v which represent respectively 
the amplitude and the group velocity of the breather, together with the space 
and time origins of the solution u v (n — riQ,t — to), one can match it to the 
asymptotic numerical solution with high precision. This enables us to identify 
correctly the emitted breathers, to determine their characteristic with a great 
accuracy, and to observe their possible disintegration into kink-antikink pairs. 

The results of such estimations are displayed in figure [7| obtained as fol- 
lows. For each driving frequency Q, we have driven the system at threshold 
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amplitude A s (£l) and then determined at the the parameters {r, v} of the emit- 
ted breathers (a given experiment at threshold driving produces repeatedly the 
same breather) . Then the amplitude of the propagating breather is plotted as a 
function of the amplitude of excitation A s (£l) for various frequencies. Similarly, 
figure [S] shows the group velocity of the breather as a function of A S (Q,). 
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Figure 7: Amplitude of the emitted breathers in terms of the amplitude of 
excitation, just above the emission threshold. 
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Figure 8: Group velocity of the emitted breathers, as a function of the amplitude 
of excitation, just above the emission threshold. 

Both figures [3 and |H1 correspond to a region of the driving frequency close to 
the gap 0.87 < Cl < 1. For lower frequencies the emitted breathers, when they 
occur, are unstable and decay into kink-antikink pairs (this will be discussed 
later). The main conclusion here is the existence of a linear relation between 
the amplitude (and the velocity v) of the generated breather and the driving 
amplitude at threshold A s . So far we have no theoretical interpretation of these 
observations. 

Another characteristic of the breather is its proper frequency Ub, related to 
the apparent pulsation period T by the usual (relativistic) relation 

-6 = ~r^= = ■ (16) 
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This frequency is plotted in figure as a function of the driving frequency f2 
and, again, we note a non trivial linear dependence. 









£.1 



Figure 9: Proper frequency of the emitted breathers, as a function of the driving 
frequency S7, just above the emission threshold. 

We now turn to the stability of the emitted breathers which, for f2 < 0.88 
decay into kink-antikink pairs. This can be qualitatively described by studying 
the binding energy W, given by 



W = 2E k - E, 



16 



VT - ^ ^ V l + r 2 J ' 
where E). is the single kink energy and Eb the breather energy 5 . 



(17) 

















0.9 





Figure 10: Binding energy of the emitted breathers, as a function of the fre- 
quency of excitation, just above the emission threshold. 

We have evaluated the binding energy W of the breathers produced in the 
simulations (by plugging in (|17|) the measured values of r and v), the result of 
which is plotted in figureliniin terms of the frequency f2 of the driving boundary. 
This figure shows that the binding energy of the breather decreases with the 
driving frequency, and goes to zero for < 0.88. Thus the breather decay in 
the frequency range [0, 0.8] is well understood from formula fTTj) . 



10 



5 Energy transmission 



5.1 Basic expression 

Nonlinear supratransmission is a process where large amount of energy flows 
through the medium. Our purpose here is to evaluate numerically this energy for 
amplitudes around the threshold value. The theoretical expression of the energy 
flow is calculated hereafter for a generic nonlinearity deriving from a potential 
energy V(u n ) (the sine-Gordon case JTJ corresponds to V(u n ) = 1 — cosu n ). 
From the energy density 

1 c 2 

H n = -u n + — (u n+1 - u n ) 2 + V(u„), n=l,---oo , (18) 
and the evolution equation follow the conservation law 

^H n + (J„+i - J n ) = 0, (19) 

with the current 

J n = -c 2 ii n (u n - u n -x). (20) 

Incorporating the potential energy resulting from the coupling of the first par- 
ticle u\ to the boundary uo, the total energy of the system reads 

oo 2 

E = Y / Hn + j{ui-u Q ) 2 . (21) 

n=l 

In our case «o(i) is the driving J5J and the chain is supposed infinite with 
u n {t) —>■ as n — * oo. 

Upon time derivation, with help of the conservation law, and using the as- 
sumed asymptotic J n — > for large n, we arrive eventually at 

^ t E = c 2 u {u Q -u 1 ). (22) 
Hence the total energy injected in the medium during time T reads 

E = c 2 dt u (t)[u {t) - ui(t)]. (23) 
Jo 

Choosing for T an integer multiple of the period of excitation makes this energy 
to vanish identically in the linear case if the driving frequency falls in the FBG. 

5.2 Numerical simulations 

In the nonlinear case, expression is computed numerically on a chain of 60 
particles with a coupling parameter c 2 = 16 and an absorbing end working over 
the last 40 particles. For a driving frequency 0.9 and amplitudes running from 
1.5 to 2.0, we obtain the figure ITT1 where the bifurcation is seen to occur for 
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Figure 11: Energy E injected in the sine-Gordon nonlinear chain for T = 140 
as a function of the driving amplitude A. 



A = 1.80, the value predicted by formula 1|10|) . This simulation has been run 
for frequencies in the range [0.2,0.9], with comparable expected results. 

Our approach stems from the existence of a breather solution of the model 
equation, allowing to determine the threshold amplitude. It is then worth won- 
dering if the process is robust against non-integrability. 

To give a partial still instructive answer to this question, we have performed 
numerical simulations of two non-integrable evolution in the same class, i.e. 
with Hamiltonian l|18Jl . 

First the Taylor truncated expansion of sine-Gordon (fifth order is kept to 
ensure a confining potential at large u n ) reads as the nonlinear Klein-Gordon 
chain: 

u n - c 2 {u n+ i - 2u n + U n -i) + u n - + ^yl4 = 0, 71=1, •••,7V. (24) 

This system is solved with the boundary driving (J2J) and the energy l|23|l is 
computed for the same parameter values as for figurc lTTl The result is displayed 
on figure IT21 




1.6 l.V 1.8 1-9 2 



Figure 12: Energy in the Klein-Gordon chain l|24|l. 

By scanning the frequency range in the gap, we have obtained that the 
process occurs down to ft = 0.7 and then disappears. We suspect that for 
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such a polynomial potential energy, at high forcing amplitude, the incoming 
wave sees an almost parabolic potential, while for low amplitude driving the 
incoming wave does feel the actual structure of the local potential. 

Another interesting non integrable Hamiltonian evolution where the local 
potential has a periodic structure is the double sine-Gordon chain 

u n - c 2 (u n+1 - 2u n + u n -i) + ~[sinu„ + sin2u„] =0, n = I,- ■ ■ , N . (25) 

Once again we have solved this system in the same situation as before and 
obtained the brutal energy flow of figure El ln this case the process holds for 
any frequency (tested from fl = 0.2) just as in the sine-Gordon case. 

300 I 1 

E «=0_9 




" OS 0.9 1 1.1 1.2 1.3 

Figure 13: Energy in the double sine-Gordon chain 

Thus we have seen that in those two cases nonlinear supratransmission does 
work which is a strong indication that it is a generic nonlinear process. An in- 
teresting question is the mechanism that generates nonlinear supratransmission, 
which is discussed now. 



6 Generating instability 

The process at the origin of nonlinear supratransmission is a nonlinear instabil- 
ity. Although it is not yet fully understood, we discuss here some of its relevant 
aspects. 

To that end we explore the properties of the sine-Gordon chain submitted 
to initial-boundary value data that precisely mimic the breather b n (t) centered 
at no — as defined in (JHJ) and which is convenient to write as 

&„(*) = 4 arctan <j> n (t) , (26) 
, . , Ac sin fit . _ 

^ = n^oTn ■ (27) 

Note that this is an exact solution in the continuous limit 1/c -C 1 and {x — 
n/c, Xc — k, ft 2 = 1 - k 2 }. 
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In order to study the behavior at the threshold, we impose the boundary 
value 

u (t) = (l + e)6o(t) , (28) 
together with compatible initial data 

u n {0) = (l + e)6„(0) , 

ii n (0) = (1 + e)b n (0) . (29) 

The boundary value at the other end of the chain (n = N) can be taken as an 
absorbing end to simulate the infinite line or as the breather value b^(t), e.g. 
to check the accuracy of the solution. 

The parameter e measures the departure from the exact solution. We are 
of course interested in what happens when e is positive. For e = we simply 
generate the approximate solution b n (t) which in the numerical simulations is 
marginaly stable (and would be indeed stable for a breather centered in no < 0). 

In the case e > 0, the numerical simulations of the sine-Gordon model imme- 
diately generate nonlinear supratransmission. From the initial-boundary value 
problem, it is natural to seek a solution as a perturbation of the breather under 
the form 

u n (t) = b n (t) +er) n (t) , (30) 
which by Q obeys at order 1 in e 

rjn - c 2 {Vn+i - 2r/„ + r) n -i) + C n i] n = eD n rj^ . (31) 

The variable coefficients C„(t) and D n (t) of this equation are given by 

C n = cos b n = 1 



(l + <^) 2 ' 

D n =i s in6 n =20 ra (i 1 ~^| )2 . (32) 

The initial-boundary value problem that goes with l|31() can be taken to 
results from l|28l) . namely 

r? (t) = b (t) , (33) 
with the related initial data. Then we observe the following fundamental facts: 

• for e < the system l|3*T|> is stable showing a long period oscillatory be- 
havior, 

• for e > the system is unstable showing an exponential growth of the 
oscillations. 

These observations are illustrated by figure ^] where we have computed the 
energy transmitted by the chain l|31|l according to formula (|23|l . as a function 
of time (the points in time are chosen as integer multiples of the period 2n/Q). 
The parameters related to this experiments are N = 100, c = 10 and the energy 
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20 40 60 

Figure 14: Energy E(T) given by (|23[1 computed for the chain 1|> . (a) expo- 
nential growth for e = 0.005, (b) linear growth for e — and (c) oscillations for 
e = -0.005. E(T) varies here on the scale [0, 2 10 4 ]. 

has been calculated for 10 points (up to t — 207r/f2). We have also displayed 
results for the linear version of system (|31|l (read with e = 0) which is marginaly 
unstable showing a linear growth of the oscillations. The exponential growth 
in the case e > is the signature of the instability which is the mechanism of 
nonlinear supratransmission. The mathematical approach of this instability for 
<{3"T|) is reported to future studies. 

7 Josephson junctions array 
7.1 The model 

A Josephson junction behaves just like a single pendulum, the rotation ampli- 
tude being replaced by the phase difference between wave functions of Cooper 
pairs jHJ- By connecting the Junctions in parallel with super-conducting wires, 
the resulting model equation is just the sine-Gordon discrete system submitted 
to damping and constant torque 0| 

% + iu n + sinii„ = J + c 2 (u n+1 - 2u n + it„_i), n=l,---,N. (34) 

Here above 7 is the constant damping along the array, J is the normalized 
intensity of the applied current and time has been normalized to the plasma 
frequency. 

We consider here only the one-dimensional geometry but it is worth mention- 
ing the coupled arrays (ladders) who revealed as a means to generate discrete 
breathers by convenient initial conditions TO] with subsequent striking experi- 
mental realizations |llj . 

The applied current can have a DC component (the so-called bias) and an 
AC driving part. The problem we are interested in is the behavior of the above 
chain initially at rest and whose first junction only is submitted to AC driving 
at a frequency f2 in the gap and intensity /. Then the model results as Q34H for 
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n > 1 and the bias J = Iq together with the following relation for n = 1: 

Hi + 7U1 + sin u\ = c 2 (u2 — iti) + /o + ^ sin (35) 

It is important to remark that the above equation can be equivalently written 
as the discrete Neumann condition 

c(u\ — uq) = sin fit, (36) 

c 

when the system (|34H is assumed to hold also for n = 1. This remark actually 
allows us to consider the continuous limit in order to determine the prediction 
of the threshold of nonlinear supratransmission. 



7.2 Bifurcation threshold prediction 

The continuous version of (|34[l and boundary value Q36|l . for the variable x = n/c 
and c — * 00 , reads as the system 

utt + ju t + sin u = I + u xx , (37) 
d x u\ x=Q = B sin fit, (38) 

where we have defined B = —I/c. This continous version corresponds to a 
long Josephson junction whose extremity x = is submitted to external micro- 
wave irradiation, the amplitude B being then related to the external magnetic 
field intensity |12j . This is a Neumann boundary condition for the sine-Gordon 
continuous equation 13711 . 

The system will then adapt the breather derivative at the boundary centered 
in — xq, namely 



sin fit sinh kxq 
fi cosh 2 kxq + (K 2 /fl 2 ) sin 2 fit' 



n 0111 j it 011111 rvu-o , on . 



where now the continuous version of the dispersion relation is fl 2 + k 2 = 1 for 
evanescent waves (due to the change of space variable, we have k = cX where A 
is defined in 0). 

Upon varying the position xq, the above expression has a maximum value for 
xq = x m given by sinh 2 nx m — 1 + k 2 /Q 2 and the related maximum amplitude 
of the derivative eventually results as the simple expression 

B s = 2(1 - fl 2 ). (40) 

This is the threshold prediction for the Neumann condition l|38|l and we are 
going now to check it on numerical simulations. 

As in section we solve the discrete system (|34|l with no constant bias 
(1 = 0), without damping in a first stage (7 = 0), for a chain of 60 particles 
with an absorbing boundary on the last 30, and with a coupling factor c 2 = 25. 
The result is displayed on figure ITsl which shows good agreement except maybe 
around the frequencies 0.33 and 0.18 as in figure^] (for the same reasons) and 
close to the phonon band, which is due to the absorbing end. 
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Figure 15: Bifurcation diagram in the (B, ft) plane. The solid line is the 
threshold expression (|40|l . The crosses indicate the lowest value of B for which 
nonlinear supratransmission occurs. 



7.3 Josephson transmission line (JTL) 

Here above the situation is that of a quasi-continuous undamped and semi- 
infinite sine-Gordon chain submitted to Neumann condition 1J38I) a the origin. 
We turn back now to the discrete case, a one-dimensional finite-length (open- 
ended) array of coupled short Josephson junctions, where the first pendulum is 
submitted to an external AC driving. 

Namely we consider the Josephson transmission line system 

Hi + 7^1 + sinui = c?(ui — U\) + Iq + /sin fit, (41) 
u n + 7u„ + sinu„ = Iq + c 2 (u n+ i - 2u n + u n -\), (42) 
u N+ i = u N . (43) 

Nonlinear supratransmission becomes here the property of the JTL to transmit 
energy under the form of kinks (or anti-kinks) as soon as the intensity I of the 
AC driving of the first junction exceeds the threshold 

I s = 2c(l - fl 2 ) (44) 

as given by the definition I = —cB and from expression (|40|l . 

A new fact here is that, for Neumann type boundary conditions as ()36|) . only 
kinks (or anti-kinks) are produced, not breathers. This is the result of the fact 
that producing a kink cost half of the energy used to produce a breather, and 
of the freedom left on the boundary value uo(t) (only the difference u\ — uq is 
prescribed) allowing full 27r-rotations. For Dirichlet boundary condition as in 
preceding sections, the precription of uo(t), e.g. by (J2J prevents full rotation. 

To illustrate this property we draw in figure [Tol a typical simulation with the 
following set of parameters 

N = 10 , 7 = 0.1 , c = 3 , I a = 0.1 , (45) 

constituting a reasonable choice for an experimental situation. At frequency 
O = 0.7, the threshold (|44|l is I s = 3.06 and for an AC driving at amplitude I = 
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3.3 the chain starts to rotate, which means that a number of kinks (elementary 
2n rotation) are generated by the AC driving. In the case of hgure ^] about 
184 kinks have been generated during 1000 units of normalized time. 



lOOO - 




o 



200 400 600 800 1 



OOO 



Figure 16: Motion of the particle 5 in a JTL of 10 junctions with parameters 

Note that we have assumed a small constant bias DC current Iq = 0.1 
in order to select kinks instead of possible anti-kinks. However the process 
does work without bias, though in some of the experiments the rotation would 
stop after some time (for a reason that we do not understand). Finally we 
mention that we have selected compatible boundary condition, namely, first an 
AC forcing that starts from zero and slowly reaches the value I within 100 units 
of time, second the following initial positions and velocities 



Then the nonlinear supratransmission does work in such JTL, despite the small 
number of oscillators, the presence of damping, the free-end boundary condition 
and the presence of constant bias. 

7.4 Energy transmission 

In order to compute the transmitted energy in the present situation, we start 
with expression (|21J) but written for a finite number N of oscillators, i.e. 



As before, this expression is differentiated and, by use of the conservation law 
<j*H3|) . the free boundary condition ujv+i = un leads to the same result, namely 



u n (0) = arcsin(io) , u n (0) = 0. 



(46) 




(47) 



71=1 



— E = c 2 u (u - ui). 



(48) 
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Now, to compute the total energy injected in the medium during time T, it 
is necessary to use integration by parts together with the Neumann condition 
(|36|l and to chose for T an integer multiple of the period of the driving. We 
eventually obtain (£ = 20 in the numerical simulations) 



„2tt 



i- 



E 



dt U\{t) cos57t . 



(49) 



This is the quantity that we have computed, for each value of the AC-driving 
/, in figure IT7I 

The bifurcation of energy translission is now quite clear on the graph of the 
energy E(T) transmitted to the chain displayed in fjgure[T7lobtained for f2 = 0.7 
and the parameters given in H45|). 
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Figure 17: Energy transmitted for 180 time units by a JTL in terms of the 
driving intensity /. 



8 Conclusion 

The ability of a nonlinear medium to transmit energy when submitted to peri- 
odic boundary driving at a frequency in a stop gap and at an amplitude beyond 
a threshold value, which we called nonlinear supratransmission, has been shown 
to have some universality (so far in the domain of nonlinear wave type of equa- 
tions). The process does not rely on integrability and is robust against damping, 
discreteness (in a reasonable range) , finiteness and different classes of boundary 
values. The mechanism at the origin of this process is a nonlinear instability 
which is still under study. 

Moreover, the nonlinear supratransmission has a simple experimental real- 
ization in the pendula chain which works surprisingly well. Some other experi- 
mental results are expected in Josephson transmission lines. 
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